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The photometry profile of the integrated Sachs-Wolfe (ISW) effect recently obtained by the Planck 
consortium by stacking patches of Cosmic Microwave Background (CMB) sky maps around a large 
number of cosmic voids, contains a cold ring at about half the void’s effective radius surrounded by a 
hot ring near the void’s boundary. The source of the temperature structure is assumed to be the ISW 
effect but the exact cause of the ringed structure is not currently well understood, particularly the 
outer hot ring. Numerical simulations have suggested that hot/cold ring structures can be produced 
by motions associated with nonlinear growths of cosmic structures whose gravitational potentials 
produce the ISW effect. We have recently developed the embedded lens theory and the Fermat 
potential formalism which can be used to model the ISW effect caused by intervening individual 
lens inhomogeneities evolving arbitrarily. This theory only requires knowledge of the void’s projected 
mass profile as a function of the passing CMB photons’ impact radius and the rate of change of that 
mass distribution at passage. We present two simple embedded void lens models with evolving mass 
densities and investigate the ISW effect caused by these lenses. Both models posses expanding mass 
shells which produce hot rings around central cold regions, consistent with the recent observations. 
By adding a small over-density at the void’s center we can produce the slight positive temperature 
excess hinted at in Planck’s photometric results. We conclude that the embedded lens theory and 
the Fermat potential formalism is well suited for modeling the ISW effect. 
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I. INTRODUCTION 

The secondary anisotropy of the Cosmic Microwave 
Background (CMB) caused by the gravitational poten¬ 
tials of inhomogeneities along the line of sight to the 
CMB photon’s last scattering at redshift z ~ 1100 is 
called the integrated Sachs-Wolfe (ISW) effect 0 . The 
ISW fluctuations caused by nonlinear growths of the den¬ 
sity perturbations at low redshifts were investigated in [2j 
using Swiss cheese models and their presence is re¬ 
ferred to as the Rees-Sciama (RS) effect. Following the 
pioneering work of 0,1 the ISW/RS effect has been 
studied extensively using analytical and numerical tech¬ 
niques such as exact general relativity (GR) modeling, 
approximate perturbation modeling, and N-body simu¬ 
lations [(i|-[23|. For example 0 studied the RS effect for 
galaxy clusters using a simple Swiss cheese model (an ex¬ 
panding homogeneous dust sphere) whereas [HI related 
some nonlinear growth effects to the divergence of the 
large scale structure’s momentum field (see also fl4l.[2ll|). 
This second order nonlinear ISW contribution is similar 
to that caused by the transverse motion of a gravita¬ 
tional lens, i.e., the Birkinshaw-Gull effect 01]. The 
numerical modeling is usually based on the merging halo 
model for large scale structure clustering [24l - |26| . Under¬ 
standing the ISW/RS effect can be considered critical to 
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modern cosmology. For example, the ISW effect probes 
the evolutionary history of cosmic structures and the dy¬ 
namics of dark energy [27| • Detecting this effect will 
directly probe the negative pressure nature of the dark 
energy and complement geometrical probes such as SNe 
la as standard candles J2|| |29| and baryonic acoustic os¬ 
cillations as standard rulers m Accurately measuring 
the ISW/RS effect is also important for quantifying the 
non-Gaussian signatures in the primordial density fluc¬ 
tuations [3lH33| . Because of the equivalence principle 
the ISW effect is manifested as a frequency independent 
temperature shift in the CMB spectrum. It has been de¬ 
tected as correlations between the CMB temperature sky 
maps and tracers of large scale structures J3J, [35| . It has 
also been detected by stacking patches of CMB sky maps 
around known large scale structures (galaxy clusters and 
cosmic voids). The ISW effect has been detected by sev¬ 
eral groups using this method 02,Q; however, the signal 
seems to be much larger than predicted by both linear 
growth theory and numerical simulations jl6H2(i 137113a ]. 
One interesting result is the strange shape of the CMB 
temperature profile of stacked cosmic voids, in partic¬ 
ular, the expected cold central region is surrounded by 
an unexpected hot ring in the outer part of the profile 
[la , I20l l3l| . This structure is hard to explain within the 
framework of the linear ISW effect in a standard ACDM 
cosmology (see Fig. 9 of HU). The stacked photometric 
profiles also show a hint of a small positive excess (at 
about 2 a) when small filter radii (< 20% of the effective 
void radius) are used. The Planck Collaboration et al. 
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na suggested small over-densities near the centers of the 
stacked cosmic voids (intrinsic to the void-finding algo¬ 
rithm ZOBOV used for the void catalog [39, .4C|) as the 
possible cause of this intriguing feature. 

We have developed the embedded lens theory in recent 
papers [41140 using the Swiss cheese models. We have 
introduced the concept of the Fermat potential (equiva¬ 
lent to the sum of the geometrical and potential part of 
the gravitational lensing time delay) for embedded lenses 
and have shown that the lowest order embedded lens the¬ 
ory can be obtained by applying a variational principle to 
the Fermat potential 1131 . By lowest order we mean small 
angle lensing caused by Newtonian perturbations of the 
background cosmology. The source of these Newtonian 
perturbations are density variations, possible large in am¬ 
plitude, but small in physical dimension compared to the 
lensing distances involved and small compared to the ra¬ 
dius of the universe. They are also assumed to be slowly 
varying. These constraints on the perturbations exclude 
the necessity of including post-Newtonian corrections 01 
in the gravity theory and result in an instantaneous New¬ 
tonian potential which satisfactorly describes transiting 
photon orbits. These restrictions are precisely the same 
as those imposed in conventional linear lensing theory. 
Our Fermat potential approach to lensing uses the pro¬ 
jected lens mass density directly without ever computing 
the Newtonian potential. 

By using the Fermat potential to formulate gravita¬ 
tional lensing the ISW effect produced by individual in¬ 
homogeneities can be obtained from a derivative with 
respect to the lens’ redshift z d of the potential part of 
the lensing time delay T p 0, 

AT „ 0T p HdT p ^r s H d f 1 dx'df(x',z d ) 

= Hd T~ = -^ 2(1 +Zd) - / —-^- 

T oz d 1 + z d c J x x' Oz d 

(1) 

where H d = H(z d ) is the Hubble parameter at the lens 
redshift z d , and the potential part of the time delay is 
defined by 

cT p (6i,z d ) = 2(1+z d )r s f ^ x 1 Zd ^ J KW ^ X ^ dx'. (2) 
Jx x 

In the above 6i is the lensing image angle, x = 9 i/6m is 
the normalized image angle, 9m the angular radius of the 
comoving lens boundary, r s is the Schwarzschild radius 
of the lensing inhomogeneity, and f[x) = M(x)/M is the 
(projected) fraction of the lens mass contained within 
the impact disk of radius 9i, and fiw{x) = 1 — (1 — 
^.2)3/2 correS p 0nc ii n g quantity for the homogeneous 

Friedmann-Lemaitre-Robertson-Walker (FLRW) sphere 
removed to form the embedded lens [461 ] . Embedded 
lensing differs from conventional lensing in that it ac¬ 
counts for the absence of the gravitational attraction of 
the removed Swiss cheese void’s mass. The fnw{x) term 
in Eq. © is a consequence of embedding on time de¬ 
lays when linear lensing is adequate. Embedding can 
produce percent size corrections to strong lensing time 


delays and image amplifications, and can produce simi¬ 
lar size changes in weak lensing shears. It also naturally 
predicts the repulsive lensing produced by cosmic voids 
[471 ]. Because the embedded lens mass is a contributor 
to the mean cosmic density its effect on passing pho¬ 
tons vanishes at impacts x > 1. For those photons that 
transit the lens the total CMB temperature fluctuation 
AT can be split in two parts, the time-delay contribu¬ 
tion ATr and the evolutionary contribution A7f, i.e., 
the two terms in Eq. ©. We have estimated the ISW 
signal caused by individual cosmic voids and clusters in 
[46( considering only the time-delay contribution A7r- 
This may have over or under-estimated the ISW signal 
depending on the sign of the evolutionary contribution, 
see Eq. (10) of 0. For example, in the case of linearly 
growing density perturbations, the evolutionary contri¬ 
bution cancels a significant part of the time-delay con¬ 
tribution. Consequently, considering only the time-delay 
contribution ATt over-estimates the total effect. On the 
other hand, for extremely nonlinear structures (e.g., for a 
deep cosmic void where the density contrast 5 is already 
approaching its lower bound —1) the evolutionary con¬ 
tribution ATs can have the same sign as the time-delay 
contribution ATt and consequently considering only the 
ATt contribution might under-estimate the total sig nal. 
Recent observations show some evidence 0, 0, 0] for 
the existence of cosmic voids much deeper than predicted 
by linear theory. There are theories of structure forma¬ 
tion that predict such deep voids 0-0]. It is important 
to study the ISW/RS effect in such strongly nonlinear 
growth periods to obtain an accurate understanding of 
what nonlinear effects are possible. In this paper, we 
present two new void lens models which include nonlin¬ 
ear growth of the lensing structures, both of which can 
produce photometric hot rings in the outer regions of the 
lens. 


II. TWO EVOLVING EMBEDDED VOID LENS 
MODELS 

Any embedded lens can be thought of as a rearrange¬ 
ment of the mass M in a comoving sphere of the homoge¬ 
neous background cosmology. We have constructed prob¬ 
ably the simplest embedded lens models in 0 to study 
the ISW effect caused by cosmic voids and galaxy clusters 
(central under and over-densities, respectively). The sim¬ 
ple void lens in [0 contains two components: an under- 
dense homogeneous dust sphere of density p = {1 — £)p 
containing a total mass (1 — £)M surrounded by a com¬ 
pensating over-dense thin shell containing the compen¬ 
sating mass . Here p is the cosmic mean density at the 
lens redshift and the parameter £ = — 6 > 0 character¬ 
izes the void’s depth: the closer £ is to 1, the deeper it is. 
Similarly simple void density profiles have been used pre¬ 
viously to model void evolution and to study weak lens¬ 
ing by cosmic voids [ 0 -[ 0 ]. For example, the inverted 
top-hat void model compensated by an over-dense thin 
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shell of finite thickness has been used in [H| to model 
cosmic voids as gravitational lenses. Very simple void 
profiles have also been used by other authors to study 
the ISW effect [II El]. To accurately compute the ISW 
effect produced by a cosmic void, linear lensing theory re¬ 
quires that we know the void’s density profile and its first 
time derivative (but no higher) at the time it lenses the 
CMB. When estimating the ISW effect caused by voids 
of such simple density profiles in j46|, we assumed two 
evolving scenarios: one where the lens was co-expanding 
with the background FLRW universe for which only the 
lensing time delay contributes to the CMB temperature 
perturbation, and a second where the density contrast 
S was evolving according to linear structure formation 
theory. The simple embedded void model developed in 
[H| predicts AT < 0 across the whole void. For that 
model there is no hot ring AT > 0 surrounding the cold 
spot as observed by Planck or predicted by numerical 
simulations In this Section we extend the model of 
[HI to include an evolutionary part, and estimate evolu¬ 
tion’s possible effect on the ISW signal. It is important 
to recall that we are measuring evolution relative to the 
expanding background Universe. 


A. Embedded Lens Model I—A Snowplow 

As with all embedded lenses the outer physical radius 
rd of the lens is related to its comoving radius \b by 
rd = R(td)xb at the time td (the time the photons en¬ 
countered the deflecting lens on their way to be observed 
at to) and where R(t d ) is the radius of the homogeneous 
universe at td . Because the mass density of the back¬ 
ground cosmology p(t d ) oc R(t d )~ 3 , the net mass M con¬ 
tained within the sphere is constant (to lowest order in 
the curvature) and equals (4/3)7rr^ p(t d )- For this partic¬ 
ular embedded void model the mass within the sphere is 
rearranged into two concentric components: a homoge¬ 
neous under-dense interior of density (1 — £)/5 and radius 
r v < I’d, and an over-dense infinitesimally thin shell at r v 
which completely compensates the under-dense interior. 
The remainder of lens mass in the outer shell r v < r < 
remains the same as the background cosmology in which 
it is embedded, i.e., is of density p[td). We assume the 
thin shell of radius r v is propagating (like a shock front) 
outward. Consequently the (normalized) radius of the 
shell y(td) = r v /r d < 1 increases with time, and the thin 
shell gains mass as it plows into p. We can also have the 
shell contract and lose mass as it fills the lower density 
inner void region, even though such motion seems rather 
un-physical. An over-dense expanding thin shell that sur¬ 
rounds an under-dense interior was proposed by several 
void formation scenarios, e.g., via explosive blast waves 
[53M55l l . The infinitesimally thin shell approximation has 
also been used to model void formation or lensing by cos¬ 
mic voids [H, H3, HH, 113] ■ The thin shell approximation 
is made for simplicity and repla cing it by a shell of fi¬ 
nite thickness is straightforward (48[. In this model any 


light ray that passes outside the compensating shell at 
y(td) is un-affected by the lens. As a consequence the 
actual size of the outer boundary of the comoving sphere 
rd is unimportant as long as r v does not overtake rd as 
the CMB photons pass through the lens. We now in¬ 
vestigate how the ISW signals produced by this model 
where y{td) = T w /r d is allowed to vary, differ from those 
of the co-expanding lens model presented in [46| where 
y(t d ) was kept constant. In what follows it is necessary to 
think of evolving lens parameters, such as the expanding 
shell’s fractional radius y, sometimes as functions of cos¬ 
mic time td at lensing and sometimes as functions of the 
lens’s redshift Zd , i.e., sometimes as y{t d ) and sometimes 
as y(zd)- The two independent parameters are related by 
IT Zd = Ro/ R(td) and obviously a quantity that increases 
with td appropriately decreases with Zd and vice-versa. 

It is straightforward to compute the projected mass 
profile of this spherical void lens 


f(x,z d )-fRw(x) = 


- £x 2 [y 2 (z d ) - z 2 ] 1/2 , 0 < x < y, 
0, y < x < 1. 


( 3 ) 

This lens profile is evolving with respect to the FLRW 
background because of the time (or equivalently the red- 
shift) dependence of the shell’s radius y. The potential 
part of the gravitational lensing time delay is obtained 
using Eqs. © and © 


cT p = 2(1 + z d )r s 



0 < x < y, (4) 


and the CMB temperature perturbation cross such a lens 
is [see Eq. ©] 


AT = 2£ r s 

T 3 c/H d 


3 yv 

XdH d 


(: y 2 -x 2 ) 1/2 {y 2 -x 2 )-^ r ,0 <x<y, 

( 5 ) 

where the redshift derivative needed in Eq. © was re¬ 
lated to the time derivative by 


d Id 

dz d (1 + z d )H d dt d ’ 


( 6 ) 


and where we have defined v = rd y{td) to be the physi¬ 
cal velocity of the shell with respect to the FLRW back¬ 
ground, i.e., the radial peculiar velocity. Both T p and 
AT vanish outside the over-dense thin shell at y{td). In 
Eq. © the two terms within the square brackets corre¬ 
spond to the time-delay and evolutionary contributions, 
respectively. They can be of the same or opposite sign 
depending on whether the shell is contracting or expand¬ 
ing. As expected this result is independent of the choice 
of r d since r s oc r[) and both x and y are oc rj 1 . Without 
any loss of generality, we take y = 1 (after Zd differenti¬ 
ation) from now on which simplifies Eq. © to 


AT 

T 


2£ r s 

3 c/H d 


(l-x 2 ) 1 / 2 




3v 

XdH d 


( 7 ) 
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and which reduces to Eq. (13) of (46| when v = 0. With 
this choice, r s is the Schwarzschild radius of the total 
mass contained within the sphere of radius r v (t d ), in¬ 
cluding the compensating shell at the time the photons 
transit the lens. Equation dTJ for AT now contains a term 
explicitly dependent on the velocity of the propagating 
shock front E GS 0 m. In Section IIIII we will use 
this equation to model the ISW cold spot across evolving 
cosmic voids possessing radially diverging (or possibly 
converging) flows. 


B. Embedded Lens Model II—Void with A 
Running Wall 


-£ 0 (y-x) 


In 




( 9 ) 


which gives an ISW effect across the lens 
AT H d T p 


T 1 + Zd 

2 r s f 
c/H d \ 


+ Q(y-x) 


2r s £ a/ 1 - x 2 /v 2 Wd) 


iv + 0 

-£0(y-x) 


In 


c/H d y H d 
1 ± yjl — X 2 4 — x 2 




ln 1 + y/1 x2 — - y/1 - x 2 /y 2 

x/y 


yj\ — x 2 jy 1 v 


y 


r d H d 


( 10 ) 


For model II we assume the lens has a physical radius 
r d = R(t d )xb and contains three mass components. The 
first is a small point mass ijM -C M at the center of the 
void (a large 77 value would produce an embedded cluster 
model, see J46|). This small (with respect to the total 
mass M of the lens) point mass was introduced to repre¬ 
sent substructures near the center of a cosmic void (e.g., a 
galaxy group or a compact object such as a black hole), 
and to possibly explain the strange small positive ex¬ 
cesses in temperature appearing at the centers of stacked 
photometric profiles [20j, [3ll |36| when filter sizes below 
^0.2 times the voids’ effective radii are used. Including 
a small central mass was also motivated by the fact that 
void-finding algorithms such as ZOBOV use galaxies as 
centers of the tessellation, and consequently the stacked 
void centers can be slightly over-dense [33, [4y] . The sec¬ 


ond component of this void model is a moving thin shell 
at fractional radius y(t d ) < 1 which contains a constant 
total mass < (1 — rj)M. The third component is a 
uniform density (1 — 77 — £)p > 0 for all x < 1 , i.e., of 
constant density contrast S = — (77 ± £). The thin shell 
can be moving out or in at the time the CMB photons 
pass, i.e., can be expanding or contracting with respect to 
the background. The mass of the shell is assumed not to 
change whether expanding or contracting, in contrast to 
model I. Without the central point mass the lens should 
behave somewhat like a cluster lens (over-density) out¬ 
side of the shell and possibly produce a hot ring there 
and behave like a void lens (under-density) near the cen¬ 
ter of the lens and possibly produce a central cold spot 
there. 

The projected mass profile of this lens is 


f(x, z d ) - /rwOe) = {v + 0(! - z 2 ) 3/2 


1/2 


(8) 

where Q(y — x) is the Heaviside step function (equals 1 
when y — x > 0 and 0 otherwise). The potential part of 
the time delay is 


cT p = 2(1 + 2 d )r s < (77 + £) 


In 


1 ± VT~. 


4 — X 2 


\/i - ■ 


where y(t d ) = v/r d and the very last term (oc v) is the 
evolutionary contribution A7/. When 77 = 0, y = 1 and 
v = 0, Eq. (fTOl) reduces to Eq. (13) of [46 ] . 

III. EXAMPLES AND DISCUSSIONS 

We now illustrate the size and shape of the ISW signals 
produced by the two lens models given in the previous 
Section [i.e., Eqs. 0 and We choose a standard 

ACDM cosmology with fi m = 0.3, Ha = 0.7 and Hubble 
constant Hq = 70kms _1 Mpc -1 , and we choose a lens 
redshift z d = 0.5. We assume the lens to be either co¬ 
expanding with the FLRW background, i.e., v = 0, or 
evolving because the thin shell is expanding or contract¬ 
ing, v ^ 0. The moving shell gains/loses mass in model 
I or has a constant mass in model II. Model I has three 
parameters: the physical radius r v = r d of the void, the 
void deepness parameter £, and the velocity v of the mov¬ 
ing shell. Model II has two additional parameters, y{t d ) 
[with 0 < y{t d ) < 1 ] which is the radius of the moving 
shell relative to the radius of the lens, and 77 , which is the 
fraction of the lens mass contained in the central compact 
object. For a model I lens we choose r d = 50 Mpc with 
v = 0, ±500km/s, or ±1000 km/s, and £ = — S = 0.5. 
For a model II lens we choose r d = 50 Mpc, y(t d ) = 0.5 
or 0.8 with v = 0, ±250 km/s, or ±500 km/s, £ = 0.9, 
and 77 = 0 (i.e., no central compact object) or 77 = 0.01 
(a compact central object containing a small fraction of 
the lens mass). 

Given the above assumed numbers we can check the 
appropriateness of having dropped all post-Newtonian 
corrections in the gravity perturbations. The projected 
perturbed void density profiles we use are smooth and 
their gravitational effects remain strictly Newtonian even 
when 6 ~ —1. For example a large void of physical 
radius ~50Mpc at redshift z = 0 constructed from a 
~2 x 10 16 Mq Swiss cheese void produces a Newtonian 
perturbation of only 4>/c 2 « 5 x 10 -5 . If the shell is 
moving with a velocity of ~1000 km/s, a post-Newtonian 
metric correction of a magnitude ~3 x 10“ 9 is gener¬ 
ated. This is the same size as (<f>/c 2 ) 2 , the next higher 
order curvature correction from GR. Such small terms 
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are never needed. A post-Newtonian gravitational cor¬ 
rection to the metric can also be generated by a pressure; 
however, with pressure estimated from the mass motions 
in our models we find a pressure to energy density ratio 
P/(pc 1 2 ) ~ 5 x 10 -5 . Therefore, just as in conventional 
lensing calculations, pressure effects are negligible. 

The results are shown in Figs. [T|[3] respectively for 
model I, model II without a central object, and model 
II with a central compact object. Besides the tempera¬ 
ture profile A T(x), we also plot the averaged signal as a 
function of the filter size R using a compensated top-hat 
filteiQ, i.e., 


(AT)(R) 


2 


' f R 

r \/2 R 

/ T(r)rdr — 

/ T(r)rdr 

J 0 

JR 


( 11 ) 


The filtered signals are shown in the right panels of Fig.Q] 
for model I, and respectively in Figs. [2] and [3] for model 
II without and with a central compact object. 

For both lens models, the total ISW signal is propor¬ 
tional to the lens mass (oc rj)), and the amplitudes of 
the signals are of the same order, r s Hd/c , which is about 
2 x 10~ 6 for rrf = 50 Mpc (a lens mass about 7.2 x 10 16 Mq 
and r s « 0.007Mpc). The relative contribution of the 
evolutionary part with respect to the time-delay part is 
characterized by v/(rdHd), i.e., the ratio of the shell’s ra¬ 
dial velocity relative to the local Hubble expansion of the 
void’s boundary, and is of order ~ 0.1 for = 50 Mpc at 
redshift z,i = 0.5 with v ~ 500km/s assuming the con¬ 
ventional ACDM cosmology. For both lens models, the 
absolute contribution of the evolutionary part is propor¬ 
tional to both the mass and velocity of the moving shell, 
i.e., ATe oc £ur s /(crd), about 2 x 10~ 7 for ss 50 Mpc, 
£ ~ 0.9, and v ~ 500km/s. Inside the moving shell, ATe 
is positive for expanding shell (diverging flow), and neg¬ 
ative for contracting shell (converging flow). This re¬ 
sult appears consistent with the Birkinshaw-Gull effect 
for CMB temperature fluctuations where a transversely 
moving concentrated lens mass redshifts CMB photons 
on the side approached bv the moving lens and blueshifts 
them on the other side (3, 0, EH fTdl. ifTlT . An expanding 
spherical shell reddens CMB photons passing just exte¬ 
rior to its boundary and blue shifts photons passing just 
interior. The converse is true if the shell is contracting. 

For the first model, the signal is dominated by the 
evolutionary part near the lens boundary where x ap¬ 
proaches 1, see Eq. 0 . Consequently, this model al¬ 
ways produces a hot ring near the lens boundary when 
v > 0. Furthermore, for expanding shells (the magenta 
curves, i.e., the top two dashed curves, in the right panel 
of Fig.0 the compensated top-hat filtered ISW cold 
spot is most significant at a filter radius R ~ 0.6 (i.e., 
r ~ 0.6rA. consistent with recent observations and sim¬ 
ulations O, [2(| 01! Hil. Predictions of the second model 


1 This compensated filter is often used by observers to reduce the 

large-scale power contamination from the primordial CMB. 


depend significantly on both the position and velocity of 
the shell. We first discuss the case without a compact 
object at the void center, i.e., rj = 0. If the shell is lo¬ 
cated close to the center of the lens (e.g., y(td) < 0.5), 
then the lens behaves like a cluster lens, see the first row 
of plots in Fig. [2] If the shell is located near the bound¬ 
ary of the lens, then the lens behaves like a void lens 
with a cold spot toward the center, see the second row of 
plots in Fig. [2] However, if the expansion velocity of the 
shell is high enough (e.g., v = 500km/s), the lens pro¬ 
duces a hot spot, AT > 0 across the entire lens, see the 
dashed magenta curve, i.e., the top curve in the bottom 
left panel of Fig. [2] For model II the motion of the shell 
only impacts regions inside the shell, i.e., the curves cor¬ 
responding to an expanding, static, and contracting shell 
fork within the shell but merge at y(td) and beyond. If 
£ is small, then motion of the shell is less important be¬ 
cause the shell contains only a small fraction of the lens 
mass and consequently the lens mass structure, i.e., fix), 
will not change significantly even if the shell is moving 
rapidly. For comparison, in Fig. [2] we also present the re¬ 
sults for the linear growth case (the cyan dotted-dashed 
curves) where we have assumed the density perturbation 
S to be evolving according to linear structure theory. For 
this case the signal is suppressed by a factor of about 4 
comparing with the static case (refer to Fig. 3 of 0 ). 

Figure[3] shows the results for void model II with a 
central compact object containing 1% of the lens mass 
(about 7.2 x IQ^Mq). In the left panel of Fig. [3] the 
spikes at the void’s center are caused by the fact that we 
have assumed a simple point mass for the compact object 
(refer to f4(| for a more complicated model). ISW signals 
averaged by applying a compensated top-hat filter show 
a small central positive excess ((AT) > 0 for filter size 
R < 0.2), see the right panel of Fig. [3] This is encour¬ 
aging given that the ISW temperature profile obtained 
by applying the aperture photometry techniques to the 
recent Planck observations and large void catalogs 0 
did contain a strange (even if with a low significance) 
positive excess near the center of the stacked images, see 
Fig. 9 of (3lj| . This result is consistent with the sugges¬ 
tion made in Planck Collaboration et al. ,3lJ that the 
positive excess may be partially caused by the fact that 
voids used for the photometry profiles contain small cen¬ 
tral over-densities. We find that a void with a moving 
shell and a central compact object can produce both the 
hot ring near the boundary of the void and a positive 
excess near the void’s center. 


IV. CONCLUSIONS 


We have recently developed the embedded lens theory 
and the Fermat potential formalism for studying the ISW 
effect caused by individual compensated inhomogeneities 
[4lM48l | . Just as with the conventional lensing poten¬ 
tial approach, the Fermat potential approach is based on 
Newtonian perturbations of the background cosmology 
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FIG. 1. ISW temperature profiles across embedded void lens models of type I are shown in the left panel. The right panel 
shows the signals averaged by applying a compensated top-hat filter, see Eq. CD- The void is at redshift za = 0.5 with physical 
radius = 50Mpc and density contrast —0.5. The respective solid, long-dashed, and dashed curves show results for shells 
which are co-expanding with (blue), expanding relative to (magenta), and contracting relative to (red) the FRLW background. 
The relative velocities shown are ±1000 km s^ 1 and ±500 km s -1 . The cyan dot-dashed curve is for the case when the shell is 
co-moving but the density contrast 5 is evolving according to linear growth theory. A central cold spot is produced for all cases 
(left panel) but linear growth produces the smallest ISW effect. Near the void boundary the time-delay contribution diminishes 
and the ISW signal is dominated by the movement of the compensating shell, which results in a hot ring for expansion. On 
the right the ISW filtered cold spot is seen to be most significant at radius R « 0.6. 


and when applied to similar lenses both approaches give 
similar results. The Fermat potential approach, how¬ 
ever, is much simpler because it uses the projected lens- 
ing mass directly and bypasses the step of computing the 
Newtonian lensing potential. Our method can be used 
to model the ISW signals extracted by stacking patches 
of CMB sky maps around known cosmic voids or galaxy 
clusters [3l|, ujfii] , or to model the CMB cold spot on the 
south hemisphere as an ISW/RS signals produced by 
large nearby cosmic voids [22], [23], [59H6lj | . The current 
modeling difficulty is the uncertainty in the mass density 
profile of voids and its dynamics. The stacked radial void 
profiles obtained from large void catalogs [13, @ built 
from galaxy redshift surveys such as the Sloan Digital 
Sky Survey [§3[ show deep interiors toward the voids’ 
centers (<5 < —0.9) for voids of radii from a few up to 
about 100 Mpc. These profiles also show compensating 
over-dense bounding ridges (see e.g., Figure 9 of 0)- 
The thickness and profile of the over-dense bounding re¬ 
gions, as well as the extent to which they compensate 
the under-dense regions are not currently very well con¬ 
strained by observations. An additional problem with 
modeling mass densities occurs because luminous matter 
as a tracer of dark matter is biased, i.e., the underly¬ 
ing void dark matter profile may be much shallower than 
estimates based on galaxy counts. Another modelling 
complication arises if suggestions that large voids tend 
to be under-compensated whereas small voids might be 


over-compensated ItI [dij are true. The universal void 
dark matter density profile of [52| based on ACDM N- 
body simulation suggests a deep void interior, S ss —0.95 
for small voids of radius ~10h -1 Mpc, and S « —0.5 
for large voids of radii ~70h -1 Mpc. These dark matter 
profiles are much deeper than those predicted by linear 
growth theory so that estimating densities based on the 
linear theory will be in error. Rather than attempting 
to model voids with complicated uncertain analytic ex¬ 
pressions or with numerically motivated density profiles 
we have applied our Fermat based calculation of the ISW 
effect, Eq. m, to two simple embedded lens models for 
the purpose of illustrating the effect evolving densities 
can have on a void’s temperature profile. These models 
possess deep voided regions and nonlinearly developing 
mass shells that produce hot rings around central cold 
spots and can thus explain recent observations found us¬ 
ing the aperture photometry technique [3ll ]3(| . The sim¬ 
plicity of the profiles (e.g., a uniform under-dense inte¬ 
rior with an infinitesimally thin over-dense shell moving 
outward or inward) combined with Eq. (JT]) allows us to 
give results in analytical form, and shed light on how 
converging/diverging flows influence the CMB observa¬ 
tions. We have assumed these void lenses to be strictly 
compensated with net mass zero with respect to the 
FLRW background. When presenting examples of the 
two simple void lens models, we have chosen voids of radii 
50 Mpc, a size not unusual according to recent void cat- 

























FIG. 2. ISW temperature profiles for embedded model II voids without a central compact object are shown in the left column 
and the corresponding filtered signals are shown on the right. The void is at redshift Zd = 0.5 with physical radius = 50 Mpc 
and density contrast 5 = —0.9. A compensating shell is located at y{td) = 0.5 or at 0.8 (top or bottom row respectively). For 
the upper row the lens behaves like a cluster lens producing a hot spot toward the center. For the lower row most of the lens 
mass is located near the void boundary and a central cold spot is produced except for the case of an expanding shell with 
v = +500kms~ 1 where AT > 0 across the entire lens. A diverging/converging flow increases/decreases the CMB temperature 
but only within the shell (i.e., the solid, dashed, and long-dashed curves merge beyond x = y(td))- A hot ring near the void 
boundary can be produced by this lens model. The filtered ISW cold spot is again most significant at a filter radius R « 0.6, 
see the bottom right panel. 


alogs 40, 6§, HH . The 50 large voids used in (Tld have a 
mean redshift of about 0.5 and an average radius of about 
100 Mpc. The assumed void density contrast, 8 = —0.5 
or —0.9 (model I and II, respectively) was largely based 
on the stacked void density profiles of [40] . While void 
redshifts and radii are provided by void catalogs, good 
estimates for compensating shell thicknesses and their 
expansion velocities are lacking. For simplicity, we have 
assumed expanding thin shells with velocities of the order 
of a few hundred kilometers per second, i.e., a magnitude 
similar to that of the peculiar velocities of galaxies. 


Expanding the analytical work given here to more re¬ 
alistic profiles and/or more complicated flows [TtI ITsl. l~iol. 
EE is straightforward but awaits observational con¬ 
straints, and may have to be done numerically. We con¬ 
clude that the embedded lens theory and the formalism 
of Fermat potential is well suited for modeling the ISW 
effect and that mass motions within voids can easily pro¬ 
duce hot and cold rings. Mass profiles with a central 
compact object can even produce a central positive tem¬ 
perature excess as observed in (3l| . 







































FIG. 3. ISW temperature and photometry profiles (left and right panels respectively) for embedded model II voids possessing 
a central compact object. The void is at redshift Zd = 0.5 with physical radius ra = 50 Mpc. The over-dense shell is located at 
y{td) = 0.8 and contains 90% of the lens mass (£ = —0.9). The central compact object contains 1% of the lens mass (about 
7.2 x 10 14 Mq). Hot rings again appear near the boundary of these large void models and the central cold spot is again most 
significant at filter radius R ~ 0.6. However, now a small central positive excess is present for filter radii R < 0.2 due to the 
point mass at the void’s center. 
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